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Abstract 

Employing the spin-flavor representation, we investigate the structures of the doubler-mixing 
symmetries and the mechanisms of their spontaneous breakdown in four types of lattice fermion 
formulation. We first revisit the U(4) x U(4) symmetries of the naive fermion with the vanishing 
bare mass m, and re-express them in terms of the spin-flavor representation. We apply the same 
method to the Wilson fermion, which possesses only the U(l) vector symmetry for general values 
of m. For a special value of m, however, there emerges an additional U(l) symmetry to be 
broken by pion condensation. We also explore two types of minimally doubled fermion, and 
discover a similar kind of symmetry enhancement and its spontaneous breakdown. 
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1 Introduction 



Lattice field theory is the most powerful theoretical framework for investigating non-perturbative 
aspects of gauge theories such as quantum chromodynamics (QCD). In this framework, while 
gauge fields are successfully incorporated without any obstacles [T], putting fermion fields on 
the lattice is not straight-forward: A naive fermion formulation gives rise to undesirable super- 
fluous fermion modes, producing sixteen flavors rather than one. The essence of this problem 
is summarized in the form of Nielsen- Ninomiya's no-go theorem [2l|3], which states that a lat- 
tice fermion action with chiral symmetry, locality and other reasonable properties inevitably 
produces one or more pairs of particles with opposite chiralities in the continuum limit. 

To overcome this doubling problem, various types of lattice fermion formulation have been 
proposed. Among them, the staggered fermion [IHB] and the Wilson fermion [T] are utilized 
frequently for QCD simulations. In the staggered fermion formulation, the number of flavors is 
reduced to four by the spin diagonalization method [7] . The staggered fermion action possesses 
the flavor- singlet vector and non-singlet axial-vector U(l) symmetries, the latter of which is 
spontaneously broken. In addition, it has several discrete symmetries as well [8]. In the case 
of the Wilson fermion, on the other hand, doublers are eliminated by adding the Wilson term 
at the cost of the chiral symmetry: The term is not invariant under the chiral transformation 
even if the bare quark mass is set to zero, so that the additive mass renormalization is required 
to realize massless or very light quarks. At the same time, the Wilson term splits sixteen 
species doublers of the naive fermion into five branches composed of one, four, six, four, and 
one fermion modes. One of the notable features of the Wilson fermion is its complicated phase 
structure [5HT^: Parity symmetry is spontaneously broken by pion condensation, depending 
on the values of the quark mass and the gauge coupling constant. 

As is seen in the above, the existence of fermion doublers on the lattice entails the symmetry 
structure different from that of continuum theory. Hence, it is important to understand the 
underlying symmetries in lattice field theory, not only for QCD simulations with these fermion 
formulations, but also for constructions of other lattice fermion formulations such as overlap 
[laiTl], domain- wall [T5Hn], staggered Wilson [T8H2T], staggered overlap [I9l|20l[22l|23] and 
minimally doubled fermions . Moreover fermion doublers could be used to economically 

simulate QCD-like theories with many flavorsjl] which might be relevant for the construction of 
techni-color theories. 

In the case of the naive fermion, it is known that the kinetic term possesses U(4) x U(4) 
symmetriej§ and it is spontaneously broken into the diagonal U(4) [TJHJ]. In the case of the 
Wilson fermion, branches other than the physical one have been taken little heed of. Recently, 



^ See, for example, refs. [33H39j 

^ For earlier discussion on the symmetry of the naive fermion, see also [40l|4T] . 
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it has been reported in the study of the Gross-Neveu model that the symmetry enhancement 
would take place at the central branch (the third branch) of Wilson-type fermions [18] . 

The aim of this paper is to shed light on the structures of underlying continuous symme- 
tries and their spontaneous breakdown in four types of lattice fermions formulation: the naive 
fermion, the Wilson fermion and two kinds of minimally doubled fermion. For this purpose, 
we rewrite lattice fermion actions in "the spin-flavor representation" |l3llll], in which the spin 
and doubler-multiplet structures of the lattice fermions become manifest. We first re-express 
the U(4) X U(4) symmetry of the naive fermion in [71112] using the spin- flavor representation. 
We then apply the same method to the Wilson fermion action, which is invariant under only 
the ordinary U(l) vector transformation for general values of the mass parameter m. We show, 
however, that an additional U(l) vector symmetry is realized by tuning m and this symmetry 
is spontaneously broken by pion condensation. Finally, we explore the Karsten-Wilczek and 
the Borigi-Creutz minimally doubled fermion and discover that a similar type of symmetry 
enhancement and its spontaneous breakdown occur. 

This paper is organized as follows. In section HI we revisit the symmetries of the naive 
lattice fermion via the spin- flavor representation. In section |3l we discuss the symmetries of the 
Wilson fermion with emphasis on the symmetry enhancement and its spontaneous breakdown. 
We also explore minimally doubled fermions in section HI Section [5] is devoted to a summary 
and discussions. Some technical details are given in appendixes. 



2 Naive fermion and Spin-flavor representation 



In this section, we first review the U(4)xU(4) symmetries of the naive fermion [3112]. Then we 
introduce the spin-flavor representation, which simplifies the identification of symmetry in the 
case of the Wilson fermion and the minimally doubled fermions. 
The action of the naive fermion is given by 

5'nf = ^ y^(^n7M^n+A - ^n+filiii^n) + m ^ . (1) 

Throughout this paper, we consider the nondimensionalized action. As is discussed in [3112], 
the kinetic term of this action has larger symmetry than the action of the continuum theory: 



^ < = exp [i (4^)r?) + e^^^T^^^ 

X 



X 



(-)f(-) 

X 



(2) 
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Here, Tx and Tx are site-dependent 4x4 matrices: 

rj) e |i4, (-ir+---+"^75, (-ir^i7M75, (-ir^^^ (3) 

rt^ e |(-i)"^+-+"n4, 75, {-iT^i,, (-1)^^17.75, , (4) 

where = ^p^^j^Up, rip^y = + and fip^^ = Yl,p^ii,v^p- Although the kinetic term is 
invariant under the transformations with arbitrary complex 9^x^ s, the link reflection positivity 
constrains 6'^'' s to be real. In other words, only if 6*^^^ s are real numbers, the transformations 
commute with the following anti-linear operation 6: 

0[^n] = i^nu-n^+l 74 , ©[^n] = 74V^n„-n4+l • (5) 

The symmetry group ([2]) is U(4) x U(4 which is broken by the chiral condensate or the : 



mass 



term down to the diagonal U(4) generated by T^"*. Therefore, there appear sixteen Nambu- 
Goldstone bosons (NG bosons) when the symmetry is spontaneously broken. The existence of 
these sixteen NG bosons is explicitly verified from the strong coupling analysis |^B5]. 

In order to understand the relation between the U(4) x U(4) symmetries and the sixteen 
doublers, it is useful to rephrase the above results using the spin-fiavor representation |13|H1]. 
We use this representation also in the subsequent sections to clarify the symmetries of the 
Wilson and minimally doubled fermions. 
Using the field Xn defined by 

Xn = iTlTlTlT^n, Xn = ^nlTlTlTlT , (6) 

we can represent the naive fermion action ([1]) as follows: 



with r/^(n) = (-1)S-<A'"'' . (7) 



^ V^^{n) Xn + m 2^ XnXn 

i,fi n 

As in the case of the staggered fermion [ISlIll], we define the field ^(N) a^ 

*(A^k/./. = E(?) ^^(^)/- ^(^k/./. = E(?) ^^(^)/- (8) 

where = or 1, Xa(A^) = X^n+a-, 1a = Ii^l2^l3^lt^ , and ia denotes the complex conjugate 
of lA- Note that we can identify the index a with a spinor index and the indices /i, /2 with 
fiavor indices of sixteen doublers. The relation between the fields ip and \l/ is given by 

^(NUhJ. = E (?) (^27v+a)/3 ' (9) 

\ Z / a,fi 



^Thc so-called "doubling symmetry" [2 is a discrete subgroup of U(4) x U(4). 

^ The spin-flavor representation becomes more complicated in the presence of gauge fields. However, because 
the structure of global symmetry remains unchanged, we here consider the free theory for simplicity. 
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\ Z / a, /i 



(10) 



In terms of "^{N), the naive fermion action ([T]) can be written as 



- ^ f (iV) (7^ ® I4 ® I4) V^^(iV) +mJ2 ^(N) (I4 ® I4 ® I4) ^(iV) 





^(iV) (75 ® 7^75"^ ® I4) Vjvl/(iV) 



(11) 



where 



5] ^'ajlj, (A)aa' (5)/,/I (C)/,/^ Vl>, 



a' 



(12) 



V^^(iV) 



^(iV + /i) -^(iV-/i) 



(13) 



^(iV + A) - 2^(A^) + ^(A^ - /i) 



(14) 
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and the superscript T denotes transposition. In ( fTTl) . the third term breaks the vector and the 
axial- vector symmetries among sixteen doublers down to the U(4) x U(4) symmetries generated 
by (14® l4(8)u(4)) and (75®7j'®u(4)). Thus, the U(4) x U(4) symmetries are part of the vector 
and the axial-vector symmetries among doublers. The reformulated action (ITT!) is quite similar 
to that of the staggered fermion [l3llll] except that eq. ( fTTl) has an additional tensor structure 
denoted by the index /2. This is manifestation of the well-known fact that the naive fermion 
is composed of four copies of staggered fermions. Although this additional tensor structure is 
trivial in the case of the naive fermion, it is nontrivial in the Wilson and minimally doubled 
fermions and is important when we discuss the symmetry of these fermions in the following 
sections. 

3 Wilson fermion 

In this section we discuss the symmetry and its breaking in the case of the Wilson fermion with 
emphasis on the effect of the Wilson term. As is well-known, the Wilson term splits sixteen 
doublers into five branches. Since numerical simulations almost exclusively use the "physical" 
branch, which contains only one massless fermion mode, the symmetry and the structure in 
other branches have not been fully investigated so far. Therefore, here we will clarify the 
continuous symmetries and their spontaneous breaking in all the branches. As a consequence, 
we will find an unexpected symmetry enhancement and its spontaneous breaking in the central 
branch. 
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3.1 Action and symmetries 

The action for the Wilson fermion [1] is given by 

S = S^f + Sw with Sw = -^"^ {-ip' 



2lpn + V'n- 



(15) 



n,fj, 



In terms of the spin-flavor representation, the Wilson term Sw is written as 

= -^E[2^(^)(l4®7>7M)*(^) + ^(A^)(l4®7>7/.)Vjvl>(iV) 

N,ti 

+^{N) (7^75 ® 75"^ ® 7m) V^^(iV)] +ArJ2 ^(N) (I4 ® I4 ® I4) ^{N) . (16) 

AT 

The first three terms in 0161] are invariant under the ordinary U(l) vector transformation, U(l)y, 
which is defined by 



(17) 
(18) 

(19) 



(20) 
(21) 

(22) 



^(A^) ^ ^'{N) = exp [1^(14 ® I4 ® I4)] ^(iV) , 
^{N) ^'{N) = t(iV) exp [-10(14 ® I4 ® I4)] , 

^ V'n = e'Vn , V'n ^ = e-'^?/^„ , 

and the site-dependent U(l) vector transformation, U(l)y, defined by 

^(A^) ^ ^'{N) = exp [10(75 ® 7J ® I4)] *(iV) , 
^{N) ^'{N) = ^{N) exp [10(75 ® 75^ ® I4)] , 

V^n ^ < = e'(-^)"^"-""*Vn , V'n ^ < = e'^-^^-^'-'-^^V'n • 

By contrast the last term in (1161) is invariant only under the U(l)y transformation. Therefore, 
the total Wilson fermion action possesses only the U(l)y symmetry for general values of m and 
r. Interestingly enough, however, the additional U(l)^ symmetry appears if m and r satisfy 
m + 4r = 0, at which the on-site terms cancel out between the mass term and the Wilson term. 
As we will show in the next subsection, this symmetry is spontaneously broken by the pion 
condensate, (ip^ip). 

3.2 Strong coupling analysis 

Now we employ the strong coupling analysis to show that there appears an NG boson associated 
with the U(l)y symmetry breaking in the presence of the pion condensate. An effective action 
for mesons in the strong coupling limit [l2l[9l[T0] can be written in general as 



S,a(M) 



Tr /(A„,^) + tr MM{n) - tr log M{n) 



(23) 
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V V V ih°'^<^ih^'i 



where Nc is the number of colors, Tr ( tr ) means a trace over color(spinor) index, and M{n) 
is a meson field. The explicit form of the function / is determined by performing a one-link 
integral of the gauge field. More explicitly we can write 

v:i = ^:;p-c+A' K!?. = -^n+A^;c, (24) 

Tr/(A„,,) = -tTf{-M{n){P+fM{n + f,){p-f), (25) 

where 4x4 matrices are specified later. In the large Nc limit, it is known that f{x) can be 
analytically evaluated as 



fix) = - 1 - In ^ ^ ^ = X + O(x') . (26) 

However, in the following part of this paper, we will approximate ][x) as /(x) = x unless 
otherwise stated because qualitative features such as an appearance of NG bosons remain 
unchanged by this approximation. 

To calculate meson masses we expand the meson field a^ 

M(n) =Ml^Y. ^''(^)rJ ' ^ e {5, P, K., v4„, T„4 , (27) 

X 

where Mq is the vacuum expectation value (VEV) of M(n), and 



I4 -p _ 75 -p _ 7a -p _ i757a p _ 7a7/3 / . a^ (^n^ 



Then the effective action at the second order of tt^ is given by 

= A^c5^[^tr(Mo-^rxMo-^ry)vr^(n)vr^(n) + ^tr(rxP,-ryP;)vr^(n)vr^(n + /i) 

n '- /i 



where 



Nc I j^y{-'p)DxY{p)Ti'{p) , (29) 



Dxy{p) = \{DxY{p)+Dyx{-p)). (30) 
Dxy{p) = ^tr(Mo-irxMo-^ry) + ^tr(rxP-ryP;)e^^^ (31) 



^In eq. (jST]) . S, P, Va, Ac and Tap stand for scalar, pseudo-scalar, vector, axial- vector and tensor respectively. 
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In the case of the Wilson fermion, M = (m + Ar)!^ = M\y1^ and = — . By taking 

Mo = crl4 + i7r75, we have 

-Mh. ± + 8(1 - r^) 

4(l-r2) ' >^W, ^^^^ 



4r2 ' IGr^a 
8r^ 



where M? = 

1 + 

As discussed in the previous subsection, at = we have an additional U(l) symmetry, 
U(l)y. Since this parameter regime resides in the parity broken phase, in which tt^ 7^ and 
< M^, U(l)y is spontaneously broken by the VEV of vr in this case. 

To compute the meson mass, we hereafter take = 1 for simplicity. Because D{p) is block- 
diagonal, we concentrate on its submatrix Dxrip) with X, F G {5*, P, A^}. Then, by setting 
p = (vr, vr, vr, vr + irrispA), we find that the S-P-Aa sector mass mspA is given by 

cosh(m5PA ) = 1 + g _ • (33) 

Note that since the transformation (122|) involves the site-dependent quantity (— i)"i+ - +"-4^ it 
is natural to expand the momentum p around (vr, vr, vr, vr). Eq. (l33l) tells us that the meson 
becomes a massless NG boson at Mw = as expected. If we use the exact form of f{x) in the 
large Nc limit, we then obtain 



cosh(™.p,) = l + ?^fMli±^, (34) 



which again shows that a massless NG boson appears at Mw = 0. 

Before closing this subsection, it is worth noting that the point M^y = corresponds to 
the central cusp in the parity broken phase, at which six fermion modes with momentum shift, 
p = (vr, vr,0,0), (vr, 0,vr, 0), (vr,0,0,vr), (0,vr,vr, 0), (0,vr,0,vr) and (0,0,vr,vr), are expected to 
appear in the continuum limit. Although we have not yet known much about the continuum 
limit for this cusp, it is expected to describe QCD with six flavors, which is still asymptotically 
free. Therefore, if an appropriate continuum limit exists, we expect the theory in the limit 
will be Lorentz-symmetric as in the "physical" branch because the Wilson fermion action itself 
possesses the hypercubic symmetry]^ In addition, the above NG boson propagator recovers 



^ Although the third term in (|16p seems to break Lorentz invariance in the continuum hmit, it is just an 
expressional artifact: The spin- flavor representation does not respect space-time symmetries of the original 
action [46] . Actually such a term cannot exist by itself because of the original translational invariance. This 
representation is not suitable for study of Lorentz symmetry although it gives good insight into other symmetries. 
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the Lorentz-covariant dispersion relation in the naive a — )■ hmit even in the strong couphng, 
which is far from the continuum hmit. We also suggest that the Wilson fermion at Mw = 
may be applicable to the simulation of the six-flavor QCD without any fine-tunings since the 
additive mass renormalization is forbidden by the additional U(l) symmetry. At this point 
{Mw = 0), while the chiral condensate a is zero, the pion condensate n is non-zero and its 
magnitude becomes maximal. 



4 Minimally doubled fermions 

Having discussed the naive and the Wilson fermion, now we move on to the analysis on mini- 
mally doubled fermions. Minimally doubled fermions are a class of lattice fermions with only 
two physical fermion modes. Continuous symmetries and their spontaneous breaking for these 
fermions have not yet been investigated enough. Therefore, in this section, we study two 
canonical examples of minimally doubled fermions, the Karsten-Wilczek and the Borigi-Creutz 
fermion, with emphasis on symmetry and its spontaneous breaking. 

4.1 Karsten-Wilczek fermion 
4.1.1 Action and symmetries 

The action for the Karsten-Wilczek fermion [23] is given by 

•S* = 5'nf + 5'kw + S^lv + "^KW ' (35) 

with 

ir ^ 

n k=l 

= ids^iJn^Ai^n, -^KW = y J] ^"^4(^A„+4 " ^n-l) ' (3^) 
n n 

where S^l^ and iS^^ counter terms of dimension three and four, respectively [28]. This 
fermion action possesses the cubic subgroup of the hypercubic symmetry. We also note it has 
CT and P invariance while each of C and T is broken [SB] . 

In the free theory, the Dirac operator in the momentum space becomes 

-Dkw(p) =^^l^l sinp^ - ir74 ^(cospfc - 1) + 1^74 sinp4 + m + 1^74, (38) 

M k 

whose zero modes appear at p = (0, 0, 0, 0) and (0, 0, 0, vr) in the absence of the mass and 
counter terms, namely m = = 0, though such a fermionic mode can be well-defined only 
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when we consider the free theory. Note that we work on the dimensionless action as shown 
in f p6|l and fl37|) : thus r, m, ^3 and all stand for dimensionless parameters. In the analogy 
of the Wilson fermion, ^3 corresponds to the mass parameter m because both are parameters 
for the dimension three terms in the action while the parameter r corresponds to the Wilson 
parameter. Thus, although the purpose of tuning ^3 is to recover the hypercubic symmetry, it 
will be tuned in a similar way to the tuning of the mass parameter m for the chiral limit in 
the Wilson fermion in Sec. |3l On the other hand, is also tuned to recover the hypercubic 
symmetry, but there is no corresponding parameter in the Wilson fermion. 

We first consider the symmetries of the Karsten-Wilczek term S'kw- In terms of the spin- 
flavor representation, S'kw is written as 

^Kw = -f 5^ [2^'(iV) (74 ® 7J7J ® 747fc) "^{N) + ^{N) (74 ® 7J7J ® 747fc) V^,^(Ar) 

N,k 

+^{N) (747^75 ® 7^75^ ® 747fc) Vfc^(iV)] + Sir ^ ^{N) (74 ® 7 J ® 74) ^(A^) • (39) 

N 

The first three terms on the right-hand side of (15^ are invariant under the ordinary axial U(l) 
symmetry U(l)yi defined by 

^{N) ^ *'(iV) = exp [i^(75 ® 7J ® 75)] *(iV) , (40) 
^{N) ^\N) = ^(iV) exp [1^(75 ® -fj ® 75)] , (41) 

lljn^^n= e^'^'^n , ^ < = ^^e'^^^ , (42) 

and the site dependent axial U(l) transformation U(1)J^ defined by 

^{N) -> *'(iV) = exp [i^(l4 ® I4 ® 75)] *(iV) , (43) 
^{N) ^'{N) = ^{N) exp [-1^(14 ® I4 ® 75)] , (44) 

in addition to the U(l)y and U(l)y, introduced in the previous section. On the other hand, 
the last term on the right-hand side of fl39|) is invariant only under the U(l)y and U(1)a 
transformations. Therefore the Karsten-Wilczek term Skw possesses only the ordinary vector 
and axial-vector U(l) symmetries, U(l)y and U(1)a- 

Next let us consider the counter terms S'^^ and S'^w- terms of spin- flavor representation, 
they are written as 

4w = i45Z^(iV)(74«)74^®74)^(iV), (46) 

N 



^Here"site dependent" and "axial" refer to original fermion fields ip and ip. 
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4w = y I] [^(^) (74 ® I4 ® I4) V4^(iV) + ^(iV) (75 ® tJts"^ ® I4) V^^(iV)] . (47) 

N 



(3) 

We notice that the dimension three counter term S'j^-yy takes the same form as the last term 



on the right-hand side of fl39|) . Therefore, the dimension three counter term is invariant under 
U(l)y and U(l)^. On the other hand, the dimension four counter term S^^l^^ is invariant under 
all of the U(4) x U(4) transformations. 

From the above discussions, the Karsten-Wilczek fermion action possesses only the U(l)y 
symmetry for general values of m, r, ^3 and ^4, and it acquires the U(1)a symmetry at m = 0. 
If (is + 3r = and m = are simultaneously satisfied, it further acquires the U(l)y and the 
U(l)^ symmetry in addition to U(l)v and U(l)^. 

4.1.2 Strong coupling analysis 

In the case of the Karsten-Wilczek fermion, we have M = ml^ + i^d-s + 3r)7j and 

|(7^ + ir74) /i = l,2,3 p- ^ / ^(7^* - i'^74) /i = l,2,3 
174(1 + 4) Ai = 4 ' \ 174(1 + ^4) Ai = 4 ■ ^ ' 

As we discussed in the previous subsection, r corresponds to the Wilson parameter while ^3 
corresponds to the mass parameter m in the analogy of Wilson fermion. Thus, in the following, 
ds will be tuned as the mass parameter m in the Wilson fermion was tuned as m -|- 4r = in 
Sec. El 

By taking Mq = al^ + ivr474, the corresponding gap equations become 

5(i±if)±(i±M..,„-^ - 0, 

3(1 - r2) - (1 + ^4)2 , , , 7r4 

J '—^ i^;r4 - (cis + 3r) - = 0. (50) 

Since the VEV of 714 would violate the Lorentz invariance, we would like to consider a solution 
with cr 7^ and 714 = 0. To obtain a solution with 714 = 0, we need to tune ^3 + 3r = for the 
dimension three counter term. In this case we need to impose 



^ = + 3(1 +r2) + (1 + ^4)2. (51) 



2 _ 2m 

As discussed in the last subsection, at m = and ^3 + 3r = 0, the action is invariant under 
an additional U(l)y and U(1)J^ symmetries in addition to U(l)y and U(l)^. Since the chiral 
condensate {ipni^n) 7^ spontaneously breaks U(l)^ and U(l)y, we expect two massless NG 
bosons at m = 0. Note that the number of fermion zero modes at 1^3 + 3r = is not two but 
does depend on the values of r and d^. At r = 1 and d^ = 0, for example, there appear six 
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zero modes at p = (0, 0, n, n/2), (0, vr, 0, vr/2), (vr, 0, 0, 7r/2), (0, n, n, —7i/2), (tt, tt, 0, — 7r/2) and 
(7r,0,7r,-7r/2). 

The inverse meson propagator matrix for S-P-Tap sector can be factorized as 



jjSPT 

where the matrix D^^ is defined by 

Daip) 

Ds{p) 



D^^ 
Z}^^ ' ' 



(52) 



SabDaip) -ifSa 

ifSfc Ds{p) 



1 1 

1 1 



2^2 4 



^{(1 + r^) _ 24, - 2r^}ck - (1 + ^4)^04 
^(1 + r^)ck + (1 + 4)^04 



(53) 
(54) 

(55) 



for X = (Ta4, S) and F = (Th^, S), and the matrix D^^ is defined by 

-abe r , 



DFvip) 

Dabip) 

Dp{p) 



i^e^'^^Se Dp{j>) 



2a2 4 



2a2 4 



^{(1 + r^) - 24a - 246}Cfc + (1 + 4)'C4 

^(1 + r2)cfc + (1 + 4)'c4 



(56) 
(57) 

(58) 



for X = {Tab, P) and Y = {Tcd,P)- Here Sa = sin(pa) and Ca = cos(pa). Since D° {n + p) = 
D^^{p), it is enough to first consider D^^{p) only and then double the degeneracy of the 
spectrum. 

By taking p = (0, 0, 0, imp), we obtain 

1 2m 



cosh(mp) 



1 + 



(59) 



This mode corresponds to a massless NG boson at m = 0. Together with the corresponding 
mode in the S sector, totally two NG bosons appear at m = as expected. 

If we introduce a small but non-zero spatial momentum, energy of the meson in the P — T 
sector becomes E"^ = rup + X(p)2 where 

2r2" 



X 



T 



mp 



1 



sinh(mp) (1 + d^) 



1 + - 



T 



ra 1 
2^ ^ 2 



{l + 3r2 + (l + rf4)'}. 



(60) 
(61) 
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To require the Lorentz covariance for the dispersion relation at m = 0, we need to choose ^4 as 

(l + dif = -r^ + VlT^. (62) 

If we use the exact form of f{x) in the large Nc limit, we then obtain 

m 

coshmp = 1 + , (63) 

'^ + *>^ . (64) 
4Vl-ff=(l+<i4)2 

which again corresponds to a massless Goldstone boson at m = 0. The introduction of small 
but finite spatial momenta results in 

= ml + X{p)\ (65) 

X = -T—^ \B-^] , (66) 



Asinhmp \ T' 
T' = 2 (a + ^) + 4(5 - C) , (67) 

B = , ^ , (68) 

2Vl - 4a2 ^ ' 
2 

C = ^ . (69) 

(1 + VI - 4a2) VI - 4(t2 ^ ^ 

This shows that we can recover the relativistic dispersion by tuning d^. As the discussion here is 
only on the strong coupling limit, the above results are not directly related to the properties of 
the continuum limit. However, the strong coupling analysis in this subsection exhibits several 
important features, such as the rotational symmetry breaking and its restoration by tuning 
counter terms, which are expected to appear also in the study of the continuum limit. We thus 
consider that our study is useful, at least qualitatively, for the investigation of the continuum 
limit of the minimally doubled fermions at ^3 + 3r = 0. 

4.2 Borigi-Creutz fermion 
4.2.1 Action and symmetries 

The action for the Borigi-Creutz fermion [23] is given by 

S = + Sbc + 'S'bc + 'S'bc ) ('''0) 

with 

-S'bc = ^^^/'n(r-7^)(V^n+A-2^/'n + ^n-/i) ' (^1) 
n, II 
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'S'bC = iCs^^nrV'n, ^bc = ^ ^ ^nr(V'n+A " ^n-/i) , (72) 
n n, fi 

r = ^E^A^' 7; = TT^r = r - 7^, (73) 

where S'g^(^ and S^^, are counter terms of dimension three and four, respectively [23130], and T 
satisfies = 1 and {r,7^} = 1. This action has the S4 subgroup of the hypercubic symmetry 
while C, P and T symmetries are broken to the combined CPT pB] . 

In the free theory, the Dirac operator in the momentum space can be expressed as 

-Dbc(p) = i b/^ sinpf, + r-f'^{cosp^, - 1) + C4T sinp^] + m + icaP. (74) 

If we consider the case with m = = 0, -Dbc(p) has one zero mode at p = (0,0,0,0), and 
another one at p = (vr/2, 7r/2, 7r/2, vr/2) if 

5^(7/. + C4P - r7;) = 2P + 4C4P - 2rP = 0, (75) 

which determines C4 = (r — 1)/2. We note there are generically two zero modes in the momentum 
space, but their positions depend on the choice of parameters. As with the Karsten-Wilczek 
fermion, C3 corresponds to m while r corresponds to the Wilson parameter in the analogy 
of the Wilson fermion. Thus C3 will be tuned as m for the Wilson fermion and ^3 for the 
Karsten-Wilczek fermion. 

First we consider the symmetries of the Borigi-Creutz term Sbc- In terms of the spin- flavor 
representation, it is written as 

= J 5^ [2^'(iV)(7.®7j7j®7.7M)^(A') + *(iV)(7.®7j7j®7.7/.)Vj*(iV) 

+^{N) (7.7m75 ® 7^75^^ ® 7.7m) V^*(iV) 
-| Yl (7m ® I4 ® I4) ^(iV) + ^(N) (7^ ® I4 ® I4) Vj^(iV) 

N,fi 

+^{N) (75 ® 7^75^ ® I4) V^^(iV) 
-ir ^(N) (7m ® 7J ® 7m) "^(N) . (76) 

N,ti 

While the first six terms on the right-hand side of fl76p are invariant under \J{l)v, U{1)a, U(l)y, 
and U(l)^, the last term is invariant only under U(l)y and U(1)J^. Therefore the Borigi-Creutz 
term Sbc possesses only the ordinary vector and axial- vector U(l) symmetry, U(l)y and U(l)^. 
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In terms of spin-flavor representation, the counter terms S^^q and S^q are written as 

N,tj. 

4c = I 5^ [^(iV)(7.®7j7>7.7M)V,vI/(iV) 

+^{N) (7,7^75 ® 7^75^ ® 7.7m) V^,^'(iV)l . (78) 



N, 11, V 



We notice that the dimension three counter term takes the same form as the last term 



^(3) 

'bc 

on the right-hand side of (!76|) . Therefore, the dimension three counter term is invariant only 
under U(l)y and U(l)^. On the other hand, the dimension four counter term, S'^q^^ is invariant 
under U(l)y, U(1)a, U(1)^, and 

From the above discussions, the Borigi-Creutz fermion action possesses only the U(l)y 
symmetry for general values of m, r, C3 and C4, and it acquires the U(1)a symmetry at m = 0. 
If C3 — 2r = is satisfled at the same time, it possesses U(l)y and U(1)J^ in addition to the 
ordinary vector and axial-vector U(l) symmetries, U(l)v and U(1)a- 

4.2.2 Strong coupling analysis 

For the Borigi-Creutz fermion we have M = ml4 -|- i(c3 — 2r)r^ and 

P+ = l{7^(l + ir) + (c4-ir)r}, p- = l{7^(l-ir) + (c4 + ir)r}. (79) 

We note C3 will be tuned as with ^3 for the Karsten-Wilczek fermion. By taking Mq = 0"l4+i-7rrr, 
the corresponding gap equations become 

2(1 + r^ + cl + Ci)a + m - „ „ = 0, (80) 

(1 + ^2 _ _ 2c4)7rr - (C3 - 2r) - — = 0, (81) 

To have vrr = 0, we need C3 — 2r = 0. In this case, a is determined by 

1 m 
C7^ a 

As discussed in the last subsection, at m = and C3 — 2r = 0, the action is invariant under the 
U(l)y and the U(1)J[ symmetries in addition to U(l)y and U(l)^. Since the chiral condensate 
{'ipni^n) 7^ breaks U(l)^ and U(l)y, we expect two massless NG bosons at m = 0. Note that 
the number of fermion zero modes at C3 — 2r = depends on r and C4. For example, at r = 1 
and C4 = —1, there are sixteen zero modes at = 7r/2 or vr for each fi. 



2 - _ . 2{l + r^ + cl + c^). (82) 
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In the strong coupling limit, there exists non-trivial mixing between mesons, which compli- 
cates the analysis on the meson mass. Fortunately, by taking C4 -|- 1 = 0, we can avoid such 
mixing. For this choice, the inverse propagator for the S oi P sector can be expressed as 



Dssip) = 
Dppip) = 




Note that to derive the above expressions we have already used C4 + 1 = 0. By solving Dxx{p + 
^x^r) = for X = S* or P with p = [p, Ex), = 1 and Ap = 0, we obtain 

2m 2 

2(cosh(Ex) - 1) = —r—^ + iP}', (85) 

a 1 + 

which has the relativistic dispersion relation at m = 0. These two modes correspond to NG 
bosons associated with the breakdown of two U(l) symmetries as expected. 



5 Summary and Discussion 

In this paper we clarify the detailed structure of continuous symmetry and its spontaneous 
breaking in four types of lattice fermion using the spin-flavor representation and the strong 
coupling analysis. 

We begin by reviewing the U(4) x U(4) symmetries of the naive fermion from the viewpoint 
of the spin-flavor representation. These symmetries can be interpreted as generalizations of 
U(l) X U(l) symmetry of the staggered fermion and are a subgroup of the U(16) x U(16) 
symmetries, which are expected to be restored in the continuum limit. We then apply the same 
method to the Wilson fermion and discover an incidental symmetry enhancement in the central 
branch of the Dirac spectrum (or equivalently, at the third cusp in the parity-broken phase). 
In this branch, on-site terms cancel out between the Wilson term and the mass term and an 
extra U(l) symmetry emerges in addition to the ordinary vector U(l) symmetry. Using the 
strong coupling analysis, we show that this extra symmetry is spontaneously broken by the pion 
condensate. We also discuss two types of minimally doubled fermions, the Karsten-Wilczek and 
the Borigi-Creutz fermions. The symmetry enhancement and its breaking also occur in these 
fermions when on-site terms are absent. Enhanced symmetries are spontaneously broken by 
the ordinary chiral condensate in this case. 

Now let us discuss possible applications of our results. First, we suggest that the central 
branch with six species of doublers in the Wilson fermion is potentially useful for the simulation 
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of many flavor QCD, especially for the simulation with six or twelve flavors. As elucidated in 
this paper, the action of the Wilson fermion without on-site terms has the enhanced U(l) 
symmetry. Since this symmetry will prohibit appearance of mass terms through quantum 
correction, we are free from additive mass renormalization in this case. In this branch, the roles 
of a and vr mesons are interchanged. The enhanced U(l) symmetry, which is broken by the pion 
condensate, is considered as a counterpart of the flavor non-singlet chiral symmetry. Although 
this U(l) is expected to enhance to U(6) x U(6) in the continuum limit of the free theory, it is 
not clear at present whether this flavor symmetry is intact also in the presence of gauge fields. 
If we could take a continuum limit of the Wilson fermion at the central cusp which restores the 
flavor symmetry, it would be valuable since it would give an alternative approach to the twelve 
flavor QCD: there is still a controversy on whether this theory is in the conformal window or 

not [33H3n]. 

Second, we comment on the implication of our results on simulations with minimally doubled 
fermions. When one uses minimally doubled fermions, one needs to fine-tune several param- 
eters to restore the Lorentz covariance, broken by the actions, in the continuum limit. To 
clarify this parameter tuning process and restoration of the Lorentz covariance, several works 
on the nonperturbative renormalization have already been done [25H5D] . However restorations 
of the Lorentz covariance in the absence of on-site terms, discussed in this paper, is not di- 
rectly relevant to the nonperturbative renormalization of minimally doubled fermions, since a 
number of fermion zero modes is actually not minimal without on-site terms. To study the 
non-perturbative renormalization of minimally doubled fermions by the strong-coupling anal- 
ysis, we need to find other tuning points at which the number of fermion zero modes is truly 
minimal. At such points, the condensate (ip'jiilj) or {^Til)) is generally nonzero. Investigations 
into such parameter regions are currently in progress. 
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A Spin-flavor representation of fermion actions 

In this appendix, we present some details on how we rewrite fermion actions in terms of the 
spin-flavor representation. 
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Let us first describe the case of naive fermions. The action of the naive fermion, expressed 
by X as in ([7]), can be recast in the following form. 



N,^l,A,B 



^Y.Xa{N)xa{N) 



+ m 

N,A 

(86) 

where and take zero or one, and the definitions of and are given in (fT3l) and 
(fT^ . In order to further rewrite the above action using the field defined in iQ and ffTOj) . the 
following formulas are useful; 

EvMm^.B + s.-,s) (I) ff ) , = E (y) , (%T I) , 

^ \ Z / a,fi \ Z / a'J{ ^ \ Z / a,fi \ Z /a'J[ j^g-jz^ 

= (7M)a,a' 5/1, /{ , 

and 



J:v,{A){6a^,,b - 6a-,, b) (f )^ (f = E (75 f 75)^ ( 



(75)a,a' (75 7;.)/i,/l • 



JA T 

2 ^'^ 



To derive fl87|) and flHHl) . it is convenient to first consider the cases when = and = 1 
separately and put them together at the end. Using these formulas one obtains the expression 
( ITTj) for the naive fermion. 

In the case of Wilson fermions, we also need to re-express the Wilson term, which can be 
represented in terms of Xa{N) as 

= ~ ^ VM)XA{Nh^^2{6A+t,,B + SA-f,,B)XB{N) - {6A+il,B - SA-i.,B)V^XB{N) 

N,ti,A,B 



+ i^A+il,B + Sa- i.,B)VlxB{N)\ +ArJ2xA{N)XA{N). 

N,A 

(89) 

To derive the expression in terms of \E', this time one needs the variants of ( IHTl) and ( IHHl) : 

j:,um..,. + >.-,.s) (I) ff ) , - E (f ) , (f , _ 

Sa,a' (7j)/i,/( 



and 

^r/^(A)(5A+A,i? - Sa-h.b) Cy) Cf) , =Y. (^5 ^^s) J ^ , , ^, 



(757/.)a,a' (75 )/!,/{ • 

17 



Applying these formulas one can re-express the Wilson term and obtain ( lT6l) . 

The actions of minimally doubled fermions can also be re-expressed by \& in a similar way 
and the resultant expressions are f l5^ and f lTB]) . 



B Strong coupling analysis for Borici-Creutz fermion 

In this appendix, some details of our analysis in the strong coupling limit for Borigi-Creutz 
fermion are given. 

Before calculating the inverse propagator, it is useful to consider the following trace formulas. 



tr(ry„p;r^,p- 

tr {TsP^TsP; 

tr(rT„,p;rT,,p- 



tr (rpP;rT,,p- 



-tr (r^„P;r^,P-) = ^ {(1 + r^)(25., -l)-cl- c] 



4 

+ + (^a/. + 56/.) 



r(l + C4) 



-tr(r^,p;rv.,p^") ^ 
-tr (rpP+rpP 



.ab/ii/ 



l+r'^ + cl + C4 



ac<Jbd- 



fi J- 1^1 / ^ ' 

(1 + r2)(l - 2^,^-2V) + (cl + r^) + (c4-r^) 

4 



(92) 
(93) 

(94) 



_j_ 

+ {5ad + Sbc - hd - 5ac) H -. {(5a/j(56c - ^bd) + ^^{^ad - 



+ 5c^l{5ad - Sbd) + 5d^(5fec - Sac)} , 

= -tr(rT„,p;r5P,") = ^^^i^(5.,-5.,j, 



-tr (rr„,p;rpP- 



4 

r(l + C4) 



(95) 
(96) 

(97) 



Since eqs. (!93l) . (196!) and (!97j) produce mixing among several sectors, which complicates our 
analysis on the meson mass, we tune C4 in order to avoid such mixings: C4 -|- 1 = 0. 
Then, an inverse propagator for the 5 or P sector is given by 



Dsip) 
Dp{p) 



2a 



1 1 + ^ 

+ ^^Z^cosp^, 



2a 



1 1 + ^ 
— ^Z^cosP/^' 



where we have already used C4 + 1 = 0. 

An inverse propagator matrix for V ot A sector is given by 



DvavM = S, 



ah 



1 1 + 
2^ ^ 4 



{2ca - ^ c„ 



1 + 



1 + 



(98) 
(99) 



{Ca + Cb), (100) 
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ab 



1 l + r^ 
so that Dv^vS^ + P) = DAaA,{p)- 
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